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Abstract. We use the generalized entropy four-current of the Miiller-Israel-Stewart 
(MIS) theories of relativistic dissipative fluids to obtain information about fluc- 
tuations around equilibrium states. This allows one to compute the non-classical 
coefficients of the entropy 4- flux in terms of the equilibrium distribution functions. 
The Green-Kubo formulae are used to compute the standard transport coefficients 
from the fluctuations of entropy due to dissipative fluxes. 

1 Introduction 

Transport coefficients, such as viscosities, diffusivities, and conductivities characterizes the dy- 
namics of fluctuations in a system. Transport phenomenon in relativistic fluids are of great 
interest in connection with the study of astrophysical conditions such as those in neutron stars, 
the study of the cosmological conditions such as those in the early universe, as well as the study 
of the relativistic nuclear collisions such as those at the Relativistic Heavy Ion Collider (RHIC) 
and at the Large Hadron Collider (LHC). 

Knowledge of transport coefficients and associated lengths and/or time scales is important in 
comparing observables with theoretical predictions. In the study of relativistic nuclear collisions 
knowledge of transport coefficients will greatly advance our current efforts and interests in the 
use of relativistic dissipative fluid dynamics in describing the observables [1121314] . Relativistic 
heavy ion collisions offer the opportunity to study interactions between hadrons over a wide 
range of net baryon density (or baryon chemical potential) and energy density (or temperature) . 
In order to study the transport properties of hadronic matter produced in such collisions one 
should extract the transport coefficients and associated length/time scales for a given model 
of interacting hadrons. Then using relativistic dissipative fluid dynamics one can study the 
sensitivity of the space-time evolution of the system and the calculated distributions of the 
hadrons to dissipative, non-equilibrium processes. In the end one needs to compare the predicted 
distributions with those observed in experiments. 



2 Second-order entropy 4-current in MIS theories 

In the MIS theories of relativistic fluid dynamics the entropy 4-current is a function of both 
the classical variables and dissipative fluxes, i.e., — s M (e, n, u M , II, q^, tt^ v ) where e is the 
energy density, n is the net charge, is the 4- velocity, II is the bulk viscous pressure, q^ is the 
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heat flux and tt^ v is the shear viscous pressure tensor. The second order or generalized entropy 
4-current in the MIS theories of relativistic dissipative fluids [5] may be written as [6l7j 

S» = s cq u» + 0q» - ^0u» (I3 n 2 - P iq v q v + faif x * uX ) 

-P(a q^n- ai q u 7r^) (1) 

where s cq (e, n) is the equilibrium entropy density, it M is the hydrodynamical 4- velocity of the 
net charge and is to be normalized such that u^u^ = 1, (3 = 1/T is the inverse temperature. The 
non-classical coefficients aj(e, n) and /3i(e, n) in Eq. |T]) are expressed in terms of thermodynamic 
integrals or can be obtained as differentiations of the equation of state p = p(a, (3) with respect 
to a and (3 (where a = (3[i with \i the chemical potential). The transition to p = p(e, n) can be 
done through the relations 

»->£■ —HI)- < 2 » 

These non-classical second order coefficients are presented in detail in Ref. [7] . From the entropy 
4-current, Eq. ([1]), the entropy density and flux are, respectively, given by 

s = tt^S" = s cq (e, n) - i/9 {foil 2 - f3 iqu q v + f3 2 -K v K vX ) , (3) 
& = A^S V = 0q» - (3((a q»n - a x q v ^ v ) , (4) 

where A^ v = g^ v — u^u v is the projection tensor with g^ v = diag(l, —1, —1, —1) the Minkowski 
metric tensor. The entropy density s = s(e, n, q^, II, 7r MI/ ) is independent of the a, while 
the entropy flux is independent of /3j. The thermodynamic coefficients /3j(e, n) > in Eq. (|3|) 
model deviations of the physical entropy density s due to dissipative contributions to S^. The 
OLiie, n) in Eq. ((4]) model contributions due to viscosity /heat couplings, which do not influence 
the physical entropy. 
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Fig. 2. The temperature dependence of the relaxation/coupling coefficients, Qi and /3i, for a pion gas. 

The non-classical coefficients cvi and fa are shown in Figs. [U [2] and [31 Shown in these figures 
are these coefficients, multiplied by pressure so that they are dimensionless, as functions of the 
ratio of pion mass to temperature z = m/T . These non-classical coefficients are not parameters 
but they are constrained by the equation of state. Thus to compute the entropy one needs 
just the knowledge of either the standard transport coefficients (in the case of a simple one 
component fluid these are the bulk and shear viscosities and the heat conductivity), or the 
corresponding relaxation times to determine the dissipative fluxes. 



3 Second moments of equilibrium fluctuations 

We shall now use the generalized entropy discussed in the previous section to obtain the trans- 
port coefficients (cf. [5]). In the local rest frame of the 3+1 formulation [5] the entropy density 
may be written as 



s(e, n, q, 27, it) = s cq {e, n) - (fan 2 + faq • q + fan : n) 



(5) 



The generalized Gibbs equation then takes the form 



9s , 9s , 9s 9s , 9s 

as = —ae + ^an + — ■ dq + T^-dll + — : an 
oe on 9q all on 



(6) 



The probability W of fluctuations of thermodynamic variables with respect to their equilibrium 
reference values is related to the second differential of the entropy by the Einstein formula [TD] 
for probability of fluctuations 

~1 



W ph exp 



-S 2 S 



(7) 
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where S is the total entropy of the system in volume V. For small fluctuations around equilib- 
rium states the change in entropy AS = S — S eq can be expanded as 

AS* (SS) cq + ^(6 2 S) cq . (8) 

Since the entropy is maximum in equilibrium, (SS) cq —0 and (<5 2 <S% q < 0. The second differ- 
ential of the generalized entropy may be derived from the Gibbs equation, Eq. ©, and when 
the resulting expression is introduced into Eq. J7]) one obtains 

W(Se, Sn, SII, Sq, 6n) « exp{ \v fe(fe) 2 + Ip^SeSn + ^(Snf 

12 I oe 2 - dean on z 

- (3 {PoSndn + ft<5q ■ Sq + f3 2 Sn : Sn) } (9) 
The second moments of a multi-variant distribution function 

W Ps exp 



--MijSxiSxj 



(10) 



where My is the matrix corresponding to the second derivatives of the entropy with respect to 
the classical and flux variables, are given by 

(SxiSxj) = (M~ 1 )ij, (11) 

where the brackets (...) denote the average over the probability distribution. In equilibrium the 
fluctuations of the classical variables e and n are uncoupled with the fluctuations of the fluxes. 
The second moments of the fluctuations of 77, q and n are given by 

(511511) = 1 



(SqMj) = ^yki , (12) 
T 

(SnijSirki) = -Tj7A ijk i , 

P2V 

with Sij the Kronecker symbol and A ijk i = (SikSji + SuSjk - 

Recall that the classical coefficients of bulk viscosity £, thermal conductivity k and shear 
viscosity n are related to the second order non-classical coefficients ft by [7] ft = Tn /(, ft = 
Tq/nT and ft = t^IIt] where the ta are the relaxation times for the respective dissipative 
fluxes. Thus the expressions Eqs. (fl"12|) relate the dissipative coefficients £, k and n to the 
fluctuations of the fluxes with respect to equilibrium. Hence these coefficients determine the 
strength of the fluctuations or vise versa the fluctuations determine the dissipative coefficients. 
The expressions, Eqs. (fT2|) for the second moments of the fluctuations allow us to compute the 
dissipative coefficients of the generalized entropy expression from the equilibrium distribution 
function and microscopic expressions of the dissipative fluxes. 

The Green-Kubo relations for bulk viscosity, thermal conductivity and shear viscosity can 
be written as 

C = ^ f ' (sn(o)sn(t))dt , 

K= T*j (Sqi(0)S qj (t))dt , (13) 
V f°° 

T]= Tj ^ l3 {Q)S^ 3 {t))dt (i^j) , 
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respectively. In Eqs. (|13[) the brackets (...) stand for equilibrium average, and no summation is 
implied over the repeated indices. The time correlation functions can be calculated from Eq. 
(fT2"]) together with Maxwell-Cattaneo relations [Hj as presented in Ref. [5] to give 

(Sn(0)dn(t)) = (T(TnV)- 1 e X p(-t/rn) , 

(S qi (0)6 qj (t)) = XT 2 (r q V)~ l 8ij exp(-t/r g ) , (14) 
(57^(0)57^ (*)) = ^lO^Ai'H exp(-t/r w ) , 

If the evolution of the fluctuations of the fluxes is described by Maxwell-Cattaneo relaxation 
equations (cf [7]), then after integration, Eqs. (T3"|) reduce to 

( = r n ^ (511(0)611(0)) (15) 

K = 'i^W)W)) 1 (16) 



V = T*~(6nj(P)8*<j(P))- ( 17 ) 



4 Hadronic Matter 



In spite of their importance, transport coefficients of hot, dense hadronic gases are still difficult 
to calculate from first principles. Progress in the study of hadronic matter transport coefficients 
is very slow, and only a calculation of transport coefficients in the variational method [11112] 
and relaxation time approximation [13] has been done. Though these models can describe some 
aspects of the properties of the hadronic matter, whether they are realistic enough or not is 
unclear. Thus, we need to investigate the transport properties of a hadron gas by using a 
microscopic model that includes realistic interactions among hadrons. This is of great interest 
in relativistic nuclear collision [14ll5j . In this work, we adopt a relativistic microscopic model, 
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Fig. 4. Shear viscosity of meson gas as a function of temperature (cf. [14] ). 



UrQMD [16] and perform molecular-dynamical simulations for a hadronic gas of mesons in a 
box of volume V. 

We focus on the hadronic scale temperature (100 MeV < T < 200 MeV) and zero baryon 
number density which are expected to be realized in the central high energy nuclear collisions. 
Transport coefficients of hadronic matter in this region should play important roles in phe- 
nomenological models. Sets of statistical ensembles are prepared for the system at different 
energy densities. Using these ensembles, the shear viscosity coefficient of a hadronic gas of 
mesons is studied as a function of temperature. 

In computing the non-classical coefficients appearing in the generalized entropy, Eq. |T]) we 
consider a system of single component gas of pions. In Figs. [I] [2] and [3] we show the relax- 
ation/coupling coefficients for a hadronic gas of pions. The equation of state is taken to be that 
of a resonance gas of non-interacting pions of mass, = 140 MeV. Knowledge of the above 
coefficients allows one to write the primary coefficients in terms of the relaxation times. Such 
relaxation times depend on the collision term in the Boltzmann transport equation, and their 
derivation is an extremely laborious task. 

We now compute the the transport coefficient of shear viscosity of a meson gas composed of 
7r, 77, to, p and 4>. We use microscopic transport model, the Ultra-relativistic Quantum Molecular 
dynamics (UrQMD) 16J, using the Green-Kubo formulas. For detailed analysis see |14) . 

Figure [4] shows the shear viscosity coefficient results from UrQMD using Kubo relations. 
As in the variational approach the coefficient grows with temperature. The UrQMD results are 
about twice those from the variational method. This might be due to the many meson resonances 
included in UrQMD while in the variational method we only have pions. Also the cross section 
parameterizations are different in the two approaches. Figure [5] shows the relaxation time for 
shear flux in a hot pion gas calculated from UrQMD by fitting the shear stress correlations. 
The dependence of the shear relaxation time on temperature is similar to the one obtained 
using variational method. The results obtained here are about a factor of two less. The reasons 
are similar to the ones given above for the shear viscosity coefficient. The use of fluctuations 
through Kubo relations has the advantage of finding not only the transport coefficients but 
also the corresponding relaxation times. In addition it is also possible to obtain the relaxation 
coefficients such as /?2- 
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Fig. 5. The relaxation time for the shear flux of meson gas as a function of temperature (cf. [14] V 

5 Conclusion 

Using the generalized entropy 4-current and the fluctuation theory together with the Green- 
Kubo formulas we presented the expressions for the standard transport coefficients in terms of 
the correlations of the dissipative fluxes and the associated relaxation times. The correlations 
of the dissipative fluxes are related to the non-classical coefficients of the entropy 4-current. 
An interesting point provided by this analysis is the reduction in the number of independent 
parameters. In the simple one component relativistic causal fluid dynamics [7] we have eight 
parameters, namely T77, r q , T n , C, k, 77, l q n, Iqw- Equations (|12jl together with two analogue 
expressions for the correlations of the heat/ viscous couplings provide five relations between 
these parameters, so that we are left with just three independent parameters, for instance 
Tu, T q , tv. In this way the macroscopic generalized entropy plus fluctuation theory give much 
more information than the standard entropy 4-flux with the same number of free parameters. 

The Green-Kubo relations are helpful in studying the evolution of the dissipative fluxes. By 
combining these relations with the generalized entropy expression provides an interesting point 
of contact between the macroscopic and the microscopic approaches. This is also evident by our 
choice of using a microscopic transport model, namely UrQMD in extracting the underlying 
transport coefficients. 
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